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For many stochastic processes there is an underlying coordinate space, V, with the 
process moving from point to point in V or on variables (such as spin configurations) 
defined with respect to V. There is a matrix of transition probabilities (whether 
between points in V or between variables defined on V) and we focus on its "slow" 
eigenvectors, those with eigenvalues closest to that of the stationary eigenvector. 
These eigenvectors are the "observables," and they can be used to recover geometrical 
features of V . 
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I. INTRODUCTION 

For many stochastic processes there is an underlying coordinate space, V. For Brownian 
motion on V the stochastic dynamics senses the structure of V through the transition rules. 
For stochastic dynamics of (say) a spin variable, a Markov process on 2 V , if the energy 
function associated with the dynamics is local and if the spin-flip rules are local there are 
again ways in which the coordinate space can be felt. 

In this article we demonstrate how in the observable representation, introduced in 0, 
HI 0, IH, the space V can be made manifest. For example, for Brownian motion on a 
(discretized) circle, the observable representation, which is a plot based on eigenvectors of 
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the matrix of transition probabilities, gives, simply, a circle. Our work may be viewed 
as part of a general effort to use spectral properties of operators to recover geometrical 
features of an underlying space. The book of Morrow and Kodaira [H[ shows how the 
slowest eigenfunctions of the Laplacian on an algebraic manifold can be used to parametrize 
that manifold. Recent work of Coifman et al. 0\ has provided specific examples of the 
use of harmonic analysis to illustrate geometric structure. It is interesting how researchers 
with different goals have converged on a common technique. Our original objective was 
the understanding of nonequilibrium statistical mechanics, of phase transitions and the 
clarification of certain fundamental questions about the second law of thermodynamics. 
Other researchers give the geometric organization of graphs as a primary goal. Nevertheless, 
in both cases spectral analysis of a stochastic process has become the method of choice. 

Although much of our previous work focused on the metric structure in the observable 
representation, the present article is concerned more with topological features (so we have 
a continuum limit in mind) and is closer to the works just cited. The exploratory aspect of 
much of what we here present, however, allows us to study issues and spaces that as far as 
we know have not been previously considered. For example, we mix manifolds of different 
dimension (line attached to a plane), we consider power sets of a set and perhaps of most 
relevance to the general representation of graphs, consider cases where degeneracy or the 
high multiplicity of eigenfunctions leads to the need for a high dimensional representation in 
order to faithfully represent an underlying space. Nevertheless, we do use metric structure 
(as in our recalling of the barycentric coordinate theorem in Sec. US, and we also establish 
a new inequality relating dynamics to distance in the observable representation (Eq. ()1U|)). 

Sec. |n]intro duces the framework and notation. Following that, in Sec. lllll we present many 
examples of spaces, V, for which the method provides faithful representations. These include 
lines connected in various ways, for example non-planar graphs for which the observable 
representation automatically requires an additional dimension for fidelity. There are also 
spaces of higher dimensions, a decorated plane and a torus. The next Section, IIV1 deals with 
stochastic process where there is an internal variable for each point of V (specifically an Ising 
model), and once again, under the right circumstances the geometry of V can be discerned. 
For some systems the observable representation demands yet higher dimension and these 
are studied in Sec. |V| Following that, Sec. |VU gives a dynamical bound on distances in the 
observable representation. Finally, in Sec. IVllI is a brief discussion of the results. 

II. THE OBSERVABLE REPRESENTATION AND ASSOCIATED NOTATION 

States of the system are given by x,y £ X, with X a space of cardinality N < oo. The 
system moves in discrete time according to a stochastic matrix R defined as follows 

R xy = "Pr(s <— y)" = Pr [State at time (t + 1) is x | State at time t is y] . (1) 

Although many of our results are more general we here restrict attention to irreducible, 
detailed balance processes. 

With detailed balance, the eigenvalues of R are real and can written 1 = Ao > |Ai| > 
| As | * * * > 0- There is a strictly positive stationary distribution, po, which is a right- 
eigenvector of R with eigenvalue 1; correspondingly there is a constant left eigenvector, 
Ao = 1 whose eigenvalue relation (A = A R) expresses conservation of probability. Be- 
cause of detailed balance, R can also be written using a symmetric matrix that we denote 
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S. Define a = diag (y/po) 0]- Then S = a~ x R(j. Call the eigenvectors of S, ip a ; then for a 
particular eigenvalue X a the left ( u A a ") and right ("p a ") eigenvectors are related to ^ a as 
follows 

Pa(x) = \/po(x)ip a (x) and A a (x) = ip a {x)/ \J po{x) (2) 

We take the ^'s to be orthonormal, inducing a different normalization for the A a 's from 
that used in [J 0, 0, HI . This normalization is used in Sec. IVI[ although the figures use the 
earlier normalization. For the examples given there is no substantial difference. 
The m-dimensional observable representation is the following collection of points 

A = {A G K m I A = (A 1 (x), A 2 (x), A m (x)) for x e X} . (3) 

To visualize this, take the left eigenvectors to be row vectors. Write the first m eigenvectors 
one atop the other. The points of the observable representation are then obtained as the 

Points of the observable representation ( j) 
A(a?i) A(x 2 ) A(x 3 ) ... A(x N ) 

\r \- "l" \r 

Ax(xi) A 1 (x 2 ) A x (x z ) ... A^xn) ( 4 ) 
A 2 {x l ) A 2 (x 2 ) A 2 (x 3 ) ... A 2 (x N ) 

A m (x{) A m (x 2 ) A m (x 3 ) . . . A m (x N ) 

The set A is produced by putting the N points, A(xk), in M. m . In jj| we focused on the 
properties of A when there was a phase transition in the system, which, as shown in [lj, 
corresponds to eigenvalue near-degeneracy for the slowest eigenvectors (those with eigenvalue 
magnitude nearest to 1). Thus if A m is very close to 1, while X m +i is not, A is a simplex 
with states in the phases at the vertices. States that asymptotically tend to the phases lie 
in the interior. The barycentric coordinates of these points with respect to the nodes of the 
simplex provide the probabilities that these points end in one or another phase. 



columns in the following array: 



Eigenvectors A\ 
(-) A 2 



A 



Prelude 

We drop the requirement that X m +i be small. This was briefly considered in where 
we illustrated that for Brownian motion on a ring, the observable representation (for m = 
2) gave a ring, while for motion on a planar region bounded by a square the observable 
representation gave a square. It turns out that this is only the tip of the iceberg: faithful 
imaging is a property of a large collection of dynamical systems, as suggested by 0] and 
demonstrated in [6( . Fig. ^ shows the much more complicated observable representation for 
Brownian motion on a fractal structure built of triangles. This image is not a picture of the 
fractal per se, but is the set of points (A%(x), A 2 (x)) for x G X, a discrete realization of the 
fractal, and A^ eigenfunctions of the associated stochastic dynamics. 

The reason this works is that dynamical proximity induces proximity in the observable 
representation. Many indications of this were given in [4j. The general idea however contains 
a subtlety: dynamical proximity in general is not symmetric (getting from x to y may 
be easier than getting from y to x), while in the observable representation distances are 
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FIG. 1: Observable representation for Brownian motion on a fractal built of triangles. This is not 
an image of the fractal coordinate space per se. Rather it is constructed from the eigenfunctions for 
the stochastic dynamics, as described in the text. On the right is a closer look. The magnification 
can be deduced from the values on the coordinate axes. 



necessarily symmetric. The inequality we will derive in Sec. IVII provides a relation to a 
symmetric dynamical construct. On the other hand when asymmetry exists it is often the 
marker of a particularly interesting class of states. We intend to discuss this in our study of 
the Sherrington-Kirkpatrick spin glass model 



III. BROWNIAN MOTION ON V: IMAGING V 



One dimension 



As indicated, for Brownian motion on a circle the m = 2 observable representation gives 
a circle. This can easily be checked analytically. The transition matrix (or its generator 
W, related to R by W = (R — 1)/At) becomes, in the continuum limit, the Laplacian. In 
the one-dimensional case, the slowest eigenfunctions are sine and cosine (with degenerate 
eigenvalues) and a circle is obtained. 

To see what happens for Brownian motion on a line one can also go to the continuum 
limit, where it is necessary to determine the boundary conditions. The fact that our discrete 
process conserves probability, meaning that there is reflection at the boundaries, implies 
the use of Neumann boundary conditions at the ends of the line. The eigenfunctions are 
essentially cosine and cosine of twice the argument. Plotted against each other, these give 
a line. The line is not straight, but the topology of the space is preserved. (See Fig. [T2lfor 
an image corresponding to the discrete walk on a line of 100 points.) 

For the next example we consider Brownian motion on a figure "Y," three lines that meet 
at a point. In Fig. [2] is the Ai-A 2 plot. Clearly the shape has been reproduced. To see how 
this has come about we show in Fig. |3] a plot of the eigenfunctions. The points (states) are 
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FIG. 2: Observable representation for Brownian motion on a figure- Y. 

numbered so that 1, 50 and 75 are the endpoints, while state-51 is in contact with state 25. 
Note that besides the functions being continuous (the apparent jump at 50 is an artifact of 
the numbering scheme) there is a condition on the derivatives at the vertex (we now discuss 
the continuum limit). Specifically, for each Ak the sum of the quantities dA^/dx must vanish 
at each vertex. This follows from integrating the Laplacian (in the general case, for which 
these considerations hold as well) across the vertex. Finally, observe that there is a discrete 
approximation to Neumann boundary conditions. 

Brownian motion on a figure eight, Fig.EJ shows that again the topology of the coordinate 
space is preserved. This remarkable interplay of only two eigenvectors can be extended, as 
illustrated for the fractal in Fig. ^ The eigenfunctions for the fractal, shown in Fig. 
produce a faithful representation of the coordinate space, although one would be hard-put 
to establish this by staring at the figure. 

A further level of complexity is obtained for a non-planar graph (in this case K 33 ). Since 
the graph is indeed non-planar, the m = 2 observable representation cannot accurately 
display the space. Fig. |U] shows two views of the 3-dimensional plot of A 3 versus A\ and A 2 . 
Another non-planar graph, the Petersen graph, is known to have many unique properties 
Hoi ]. For the observable representation it presents features not encountered in the examples 
given so far. As usual, we associate a matrix of transition probabilities with this 10-node 
graph by taking the adjacency matrix to be equal to (the non-diagonal part of) the generator, 
W, of the stochastic dynamics. W's diagonal is fixed by W xy = 0. (The generator is 
related to "i?" by R = 1 + WAt, with At small enough for i?'s matrix elements to be 
non-negative.) 

To illustrate the new feature present in the Petersen graph, in the left-hand image of 
Fig. [7| we show the result of straightforward diagonalization and plotting. The lines drawn 
are not part of the observable representation, but only indicate which nodes are connected. 
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FIG. 3: Eigenfunctions on which Fig. [21 is based. On the left is the slowest eigenvector (eigenvalue 
of W = -0.000985). On the right is the next slowest (eigenvalue -0.00141). The fact that the 
trigonometric functions need not be the same on different branches leads to the crowding of points 
near the ends of the lines in Fig. |21 
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FIG. 4: Observable representation for Brownian motion on a figure eight . 



The reason for this hodgepodge is that there are five degenerate eigenvalues, and which 
axes are chosen in this space depends only on vagaries of the diagonalization technique. In 
the same figure, on the right, we selected within this 5- dimensional space eigenfunctions 
that had particular properties with respect to the symmetries of the graph. The permuta- 
tion (2,3,4,5,1) (7,8,9,10,6) commutes with W and has eigenvalues that are the 5 th roots of 
unity. The real and imaginary parts of the complex simultaneous eigenvectors of W and the 
permutation are then eigenvectors of W. These eigenvectors were used for the right-hand 
illustration in Fig. [7| This is one of the standard representations of the Petersen graph, and 
in particular makes manifest the 5-fold symmetry. 
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FIG. 5: Eigenfunctions on which Fig. ^ is based. 



B. More than one dimension 



Our examples so far have dealt with one- dimensional coordinate spaces. But the method 
has no such restriction. As shown in the observable representation for Brownian motion 
on a square is — simply — a square. In the next example we mix dimension. Suppose that 
there is a string attached to the square at one of its corners and the particle can wander off 
onto that string. Fig. |Hl shows the observable representation for this stochastic process. We 
remark on a feature that can be seen in this figure as well as in Fig. IHland in Fig. 11 (on the 
right) in 4]. There is a crowding of points near the edges. This is because one is plotting 
different trigonometric function against one another, and although no change in topology 
occurs, distances are not preserved. 

A 2-dimensional space requiring embedding in three dimensions is the torus (i.e., we 
consider Brownian motion on the torus). This further illustrates a limitation of our method: 
when degeneracy appears, particular mixtures of the eigenvectors may cause the observable 
representation to collapse or project degrees of freedom of the coordinate space. (This 
problem was expressed a bit differently for the Petersen graph.) In Fig. El we show a torus 
discretized at different levels for its two generating circles. For this situation, Ai = A2 
and A 3 = A 4 , but they do not equal each other. The figure shown is obtained by selecting a 
random 3-space within the four dimensions spanned by (Ax, A 2 , A 3 , A4) (and then projecting 
onto a plane for publication). However, if one takes Ai and A 2 and any direction 

in the plane (A3, A4) one gets the image shown in Fig. [TU] (provided the discretization is 
different for the two radii: otherwise there is a four-fold degeneracy, effectively randomizing 
the axes). 

Many other figures can also be constructed, for example the edge of one plane meeting a 
line through another. But these do not present problems and we omit them. 



8 



1s 




FIG. 6: Two views of the observable representation for Brownian motion on a non-planar graph. 
The dimension "m" of the observable representation must be set to 3 for fidelity. All vertices have 
3 edges. The line crossings in both figures are artifacts of the two-dimensional projections, as can 
be ascertained by comparing the figures. 




FIG. 7: Observable representation for the Petersen graph. On the left is the result of a straightfor- 
ward diagonalization using matlab. On the right the degenerate eigenvectors have been selected 
so as respect a five-fold symmetry of the graph. 
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FIG. 9: Brownian motion on a torus. The axes are three random orthogonal vectors in the 4- 
dimensional observable representation. 
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FIG. 10: Same as the previous figure, but the axes are (Ai, A2, A3). 

IV. BROWNIAN MOTION ON 2 V : IMAGING V 

Stochastic dynamics often describes the behavior of a collection of objects located in a 
space and having internal degrees of freedom whose interaction depends on the proximity of 
those objects within the space. An example is the Ising model with Glauber or Kawasaki 
dynamics. For this circumstance the path from dynamics to coordinate space is less direct. 
Nevertheless, since the dynamics does contain spatial information, appropriate probes can 
lead to imaging of the coordinate space using the observable representation. 

In Fig. ^2 we show the observable representation for a ring of particles having ferromag- 
netic nearest neighbor interactions. We use Kawasaki dynamics at low temperature |ll| . 
The number of states is of course large, and because there is a conceptual distance between 
the space "V" and the dynamics on 2 V a bit of interpretation is in order. For this dynamics, 
the stationary state, po is quite non-uniform and for those states having the largest po-weight 
we have used a larger symbol. We took 12 spins at low temperature with the condition that 
exactly 3 have the value +1. The 12 heavy points in the figure correspond the 12 possible 
clusters of 3 spins. They forms a circle because the stochastic process by which one can pass 
from one to the other of this states is ultimately a Brownian motion on a circle. That same 
symmetry exists for other possible spin configurations, 2 up, 1 down and all 3 separate. 

A variation on this theme occurs for Kawasaki dynamics on a line. To orient the reader, 
we first show what the observable representation gives for Brownian motion on a line. We 
earlier indicated how this is obtained analytically, but it is helpful to see the image. This is 
Fig. The observable representation for Kawasaki dynamics on a line is given in Fig. [T3l 
The linear nature of the model emerges through the identification of points in the plot. The 
two points with maximal weight are on the lower left and lower right. These are clusters of 
three attached to the ends of the line. Because of the boundary condition they have only 
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FIG. 11: Observable representation for the stochastic dynamics of nearest-neighbor Kawasaki 
dynamics on a ring of Ising (ferromagnetic) spins. Heavier symbols are used for states having 
greater po-weight. 
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FIG. 12: Observable representation for Brownian motion on a line. The analytic form was given 
earlier; this image is presented for comparison with Fig. 1131 



one broken bond and are favored over all other states. Clusters of 3 in the middle of the 
line require two broken bonds and are the next most likely states. These form the arc in 
the upper portion of the figure. This is because the transition from one to another is quite 
similar (if requiring more steps) to the simple walk on a line. Finally there are the two points 
near the bottom, between the line-end 3-clusters. They represent states with a cluster of 
two (up spins) on one end of the line and a cluster of one on the other. Because of the 
boundary condition, this costs only two broken bonds, giving these states the same weight 
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FIG. 13: Observable representation for the stochastic dynamics of nearest-neighbor Kawasaki 
dynamics on a line of Ising (ferromagnetic) spins. 

as internal clusters of three. These are col points, in the sense that the quickest way to go 
from having one (maximal) cluster on the left, to one on the right, is via spin flips that take 
as intermediate values these two states. In line with the barycentric coordinate theorem 
of IH, their positions give the probabilities of reaching the minima. The shape of the arc 
associated with the random walk on the line is that same as what one gets from Brownian 
motion on a line, which is curved for reasons mentioned above. See Fig. ED There is a sharp 
intrinsic distinction between the intermediate points on the bottom and those of the arc (so 
this behaves quite unlike the circle shown earlier). As the number of spins changes, the two 
intermediate points on the bottom do not change. On the other hand, increasing (say) the 
number of spins will increase the number of heavy points on the arc. 

V. PROLIFERATING DIMENSIONS 

When a graph has many loose ends, the dimension required for a faithful representation 
may become large. We will look at two examples, the star and the Bethe lattice (complete 
binary tree). The difficulty is that the extremities of the graph have eigenf unctions that are 
entirely supported there, so that a few slow eigenvectors will not separate these points. The 
Bethe lattice is often described as being of infinite dimension, and this characterization is 
supported by the observable representation as well. 

A star consists of a central point surrounded by n other points, all of which communicate 
only with the central point. The generator of the stochastic dynamics is an (n + 1) by 
(n + 1) matrix whose only off-diagonal elements are ones in the first row and column. (The 
diagonal is adjusted so that W xy = 0.) The eigenvectors to be used in the observable 
representation are (n — l)-fold degenerate. Calling them they have the property that 
their zero component is zero and the sum of components 1 through n is zero. There are two 
other eigenvectors. One is constant on all points. The other is constant on the non-central 
points, but takes a different value on the central point. As a consequence the points of the 
observable representation satisfy a single linear constraint but are otherwise independent. 
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That constraint is that their sum is the zero vector, which is true because the components 
of the individual eigenvectors sum to zero. 

We now present a specific construction showing that the observable representation gives 
a simplex in MJ 1 . Let v = n — 1. The unnormalized but orthogonal eigenvectors (reading 
horizontally, as in Eq. (J3J) can be taken to be 
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(5) 



There are v = n — 1 rows in this array and n + 1 columns. For the (n + 1) points of 
the observable representation one reads down. The first column gives the (interior) point 
corresponding to the center of the star. We must show that the other n points (i.e., columns 
2 through (n + 1), considered as vectors in M n ) have non-zero volume. It is not difficult 
to show that the v x v matrix consisting of columns 2 through n has determinant unity 
12l |. This proves that these v points have non-zero volume with respect to the origin of 



coordinates. To see that they have non-zero volume with respect to another point in W 1 
one should subtract that point (a column vector) from every one of them, and evaluate the 
resulting determinant. Now from the determinant of columns 2 through n one can subtract 
any sum of multiples of columns without changing its value. We therefore subtract the sum 
of all the v columns, which one notes is precisely the final column in the array. This proves 
our assertion 

If one produces the observable representation for any smaller dimension one still gets a 
simplex but it is not a full image of the graph: points will be coincident. 

Similar results are obtained for the Bethe lattice. Here too there are eigenfunctions with 
support entirely on the fringes, so that omitting them yields a figure in which distinct points 
of the graph coincide in the observable representation. In Fig. El we show the result of a 
3-dimensional plot of the 15 points of a tree with 3 levels of descendants. 

It is interesting that even if one fills in the lines of the graph, a faithful image is still not 
obtained. Edges of the graph were replaced by a sequence of points on which deterministic 
or Brownian motion took place. This made no substantial difference. On the contrary, this 
raises the question with respect to our Brownian motion figures: what is the minimal graph 
that would recover the same topological features. 



VI. RATIONALE 

In previous work we gave many relations of use in interpreting distance in the observable 
representation 0, 3- The following (new) relation, Eq. (jl(J|) . provides a general inequality. 



It is not particularly sharp, but applies throughout the set A, not merely at its extremities. 
Our proof only applies when R satisfies detailed balance. 

Define p x (u,t) to be the probability distribution at time t of a system whose time-0 
position was x G X. Thus p x (u,t) = R l ux . A reasonable distance between two points x 
and y could be taken to be the variation distance between these two distributions. We look 
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FIG. 14: Observable representation for a Bethe Lattice of 15 points. Note that points at the 
extremities of the graphs are not distinguished, reflecting the fact that the first three eigenvectors 
take the same values on these points. As for the star graph, increasing numbers of dimensions are 
required for a faithful representation. 



instead at something closely related: 

D(x,y;t) = J2 



Px(u,t) -Py(U,t) 



(6) 



with po the stationary state. By straightforward manipulation we rewrite R in terms of S 
and a (see Sec. |H} and use a spectral expansion for 5*. First note that 



Px(u,t) _ 



u 



X 



x 



(7) 



where use has been made of Eq. (J2J) and the definitions preceding it. We next use Eq. (J7J) 
combined with the fact that the L 1 norm exceeds the L 2 norm: 



D(x,y;t) > 



P X (u,t) -Py(u,t) 



Po{u) 



\ u 



(8) 



This expression is of the form E u S a V'a( M )V ; o'(' u ) c a c a'] 1 ^ 2 ) which because of the or- 
thonormality of the ip's is [J2 a l c a| 2 ]^ 2 - ^ follows that 



D(x,y;t) > /^|A a |*K(x)-A a (2/)f 



(9) 



The sum on the right can be truncated at the m th eigenvalue, preserving the direction of 
the inequality. Since the magnitudes of the eigenvalues decrease monotonically (with index) 
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we obtain, finally, 



m 



D(x,y;t)/|A m |* > 



\ 



^2\M*) - Mv)\ 2 ■ 



(10) 



a=l 



The quantity on the right is the distance in the observable representation. Thus, two points 
x and y which are adjacent dynamically — in the sense that the Markov chain soon nearly- 
forgets which point it came from — will be spatially adjacent in the observable representation. 



Much of our previous work on the observable representation focused on metric properties. 
We used an assumed near degeneracy among the largest (slowest) eigenvalues to establish 
that points in separate phases formed tight clusters and were far from one another (in the 
observable representation). In this article we drop the assumption of near-degeneracy. The 
information obtained thereby is therefore mainly topological. For this reason the bound 
obtained in Sec. |VJ may be expected to be useful. 

A second assumption that we drop is that after the eigenvalues used in the observable 
representation there is a significant decrease in the magnitude of the subsequent eigenvalues. 
From Eq. (fTUj) it is clear that one does not want to go to m such that |A m | is small, but not 
having |A m +i| be small presents no difficulties (whereas in jj| it was important that there 
was a range of t values for which A^ ~ 1, while A^ +1 was small). 

The message to be drawn from this work is that the topological features of an underlying 
space for stochastic dynamics are encoded in the eigenvectors. The observable representation 
is a way to read that code. In some of our earlier work la] that "code" was deciphered by 
creating clusters based on a dynamical metric (see also |3[). Here though we go directly to 
the eigenvectors. This same technique has been used by Coifman et al. and anticipated 
by Kodaira [5]. Our examples also deal with spaces that are neither manifolds nor graphs 
and also treat the recovery of geometry from a stochastic process in what is essentially a 
field defined on the underlying space. Finally, we remark that deducing spatial properties 
from the spectral properties of an operator defined on the space has a long history, with 
a well known example being Kac's "Can one hear the shape of a drum?" [ljjjj]. There 
however it is asymptotic properties of the eigenvalues that provide the information. Here it 
is the eigenfunctions of the slowest eigenvalues, functions that for physical reasons we term 
observables. 



We thank Bertrand Duplantier, Thomas Gilbert, Peter Greiner, Annick Lesne, and Jean- 
Marc Luck for helpful discussions. This work was supported in part by NSF, NSA and ARO 
grants. 



VII. DISCUSSION 



Acknowledgments 



[1] B. Gaveau and L. S. Schulman, Theory of non- equilibrium first order phase transitions for 
stochastic dynamics, J. Math. Phys. 39, 1517-1533 (1998). 



16 



[2] L. S. Schulman and B. Gaveau, Coarse grains: the emergence of space and order, Found. Phys. 
31, 713-731 (2001). 

[3] L. S. Schulman and B. Gaveau, Dynamical distance: coarse grains, pattern recognition, and 
network analysis, Bull. Sci. math. 129, 631-642 (2005). 

[4] B. Gaveau and L. S. Schulman, Multiple phases in stochastic dynamics: geometry and proba- 
bilities, Phys. Rev. E 73, 036124 (2006). 

[5] J. Morrow and K. Kodaira, Complex Manifolds (Oxford University Press/American Mathe- 
matical Society, Oxford, 2006), original, Holt, Rinehart and Winston, New York, 1971. 

[6] R. R. Coifman, S. Lafon, A. B. Lee, M. Maggioni, B. Nadler, F. Warner, and S. W. Zucker, 
Geometric diffusions as a tool for harmonic analysis and structure definition of data: Diffusion 
maps, Proc. Natl. Acad. Sci. 102, 7426-7431 (2005). 

[7] In using the operator "diag" we follow matlab's notation. Acting on a vector, diag creates 
a diagonal matrix with that vector on the diagonal, zeros elsewhere. Acting on a matrix, it 
produces a vector whose components are the diagonal of the matrix. 

[8] B. Gaveau and L. S. Schulman, Mean Field Spin Glass in the Observable Representation 
(2006), in preparation. 

[9] D. A. Holton and J. Sheehan, The Petersen Graph, vol. 7 of Australian Mathematical Society 
Lecture Series (Cambridge U. Press, Cambridge, 1993). 

[10] M. Aigner, Combinatorial Theory (Springer- Verlag, Berlin, Heidelberg, New York, 1979). 

[11] In Kawasaki dynamics the random process allows (only) nearest-neighbor double spin flips for 
oppositely directed spins, hence preservation of the total magnetization. 

[12] One can show that the determinant is unity by induction (expanding about the right hand 
column of the matrix and doing successive expansions). Alternatively, recall that the rows in 
Eq. (J5J are orthogonal. Hence the determinant is the product of the Euclidean norms of the 
rows. Those are Vl + 1, y 7 ! + 1/2, \J\ + 1/3, . . . , except for the last row. Its value is \J\jv 



[13] A second proof of the same assertion is as follows. First, none of columns 2, . . . , n+ 1 in Eq. (J5J) 
is within the convex hull of the others. Regarding columns 3 through n + 1, each of them has 
a negative coefficient in a location different from that of the others. As to column 2, it is not 
within the convex hull of the others, because any convex combination of the remaining columns 
must have a negative coefficient in the coordinate corresponding to the first column used with 
positive weight. The same argument also shows that any smaller subset of the columns are 
lower dimensional faces of the larger simplex. 

[14] B. Gaveau and L. S. Schulman, Creation, dissipation and recycling of resources in non- 
equilibrium systems, J. Stat. Phys. 110, 1317-1367 (2003). 

[15] M. Kac, "Can one hear the shape of a drum?", Amer. Math. Monthly 73, 1-23 (1966). 




